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Abstract 

Lotka Volterra model and its modified forms have long become a major 
area of interest for periodic motions in nonlinear systems with compet- 
itive species. The model given by Volterra shows that its periodicity is 
dependent on initial condition. This characteristics allows us to calculate 
the effect of periodic seasonal changes on population densities of different 
species. 

1 Introduction 

In the third decade of the last century Lotka £Q and Volterra [2] formulated a 
coupled set of equations to describe an auto catalytic model and the statistics of 
fish catches in Adriatic. Since then Lotka Volterra model has become a central 
area of interest for periodic oscillations in nonlinear systems with competitive 
elements |3J< It has found its applications in population biology 0], ecology 
mathematical biology [7]. Nevertheless It has often been criticized for 
being biologically unrealistic and mathematically unstable. Plenty of modifi- 
cations were made to make it more realistic. Yet in last two decades area of 
application of Lotka-Volterra model expanded. It has made its way through 
newly explored applications, from membrane dynamics of competing neurons 
H3, neural networks [Sj , metabolic algorithm JOj to network-electronics |11II12| 
and stochastic dynamics ^21- With this re-emergence of Lotka-Volterra dynam- 
ics, we study some unexplored aspect of the original model. 
Lotka-Volterra dynamics describes a predator-prey system where the prey thrives 
on the naturally available resources while predator thrives only on its interac- 
tion with the prey and would be extinct in a non interacting system. It is 
easy to realize that in such a system the prey population would thrive of the 
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predators decrease while an increase in the predator number would cause the 
prey population to diminish. The predator cannot flourish for ever because 
a decreasing number of prey would lead to fewer interaction and would drive 
the predator towards extinction. This produces the limit cycle like oscillations. 
The Lotka Volterra limit cycle is however very different from the usual limit 
cycles in that it depends on the initial conditions. This makes Lotka Volterra 
dynamics interesting. The unbounded growth of the prey population in the 
absence of the predator depends on a constant source of nutrients. In any prac- 
tical situation, there will be some periodic fluctuations in the nutrients because 
of natural causes. If it is a wolf-rabbit situation, then the rabbit's supply of 
grass will have a seasonal variation. Interestingly enough, the periodic variation 
in the growth rate of the prey has not been studied in literature. There has 
been study of random variation |16j but not periodic ones. In this work, we 
study the effect of periodic variations. In carrying out this analysis, we noticed 
that the standard techniques that are used in dealing with nonlinear oscillators 
(Poincare Lindstedt method 14 , equivalent linearization ^j] have not been car- 
ried out for the Lotka volterra model. Accordingly in section 2. We apply the 
Lindstedt Poincare technique on the Lotka Volterra model and come up with 
an initial condition dependent oscillation period. Our analytic and numerical 
results agree. In section 3, we use the analytic tools of section 2 to study the 
dynamics of the modulated system and also carry out a numerical analysis to 
test our predictions. Since thcunmodulatcd system shows an initial condition 
dependent limit cycle, the final state of the driven system is found to be initial 
condition dependent. In section 4 , we apply another specified technique which 
work for high frequency perturbation and show how the calculation numerics 
compare favourably. A brief conclusion is presented in section 5. 



2 Depencence of Initial Condition 

In this section, we focus on the limit cycle of the Lotka- Volterra model. As is 
well known, this limit cycle is unusual in the sense that it depends on the initial 
conditions. We show how the Poincare-linstedt technique can be applied to the 
model to obtain an initial condition dependent frequency for the limit cycle. 
This calculation is valid for small amplitudes and our numerical calculation 
supports the validity of this result. 

We take the simplest form of two species Lotka- Volterra system; 

x = x — xy , 

(2.1) 

y = -y + xy 

According to a linear stability analysis, the fixed point at (1,1) is a center. We 
first transform Ea H2.1|) to the variables Xx = x — 1 and yx = y — 1 to write, 

xi = —3/1 — xiyi 
y\ =x 1 + xiyi 
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We note immediately that if we drop the nonlinear terms, then we have a simple 
harmonic oscillator of frequency unity. We are now in a position to carry out 
a Poincare-linstedt analysis. We imagine the existence of a parameter A multi- 
plying the non-linear terms in Ea i|2.3[l and introducing the frequency u> of the 
full dynamics, rewrite Ea i|2.2[l in the form; 

K:iMi *)(sM (-";)<-> 

For A<1 we expand 

oj =1 + Aa>i + X 2 uj2 + ... 

xi = %w + Axu + \ 2 x\2 + ■ ■ ■ (2.4) 

y\ = yio + Ayn + X 2 yi 2 + ■■■ 

The right hand side of Eg (|2.4[l has term of order A and higher. At order unity, 
the solution is xw = ReAe lult ,y±o = Re4-e tuJt , where A is the amplitude of 
motion. At 0(A), we have 

'( In H'l 1 + -to) <") 

The solvability condition for an inhomogeneous second order differential equa- 
tion now leads to uj\ = 0. The solution for x\\ and yn is found to be 

A 2 



x n =Re-^-(l + 2i)e 2t " t 
yn = Re^(l-2i)e 2i ^ 



(2.6) 



AtO (A 2 ) 

L ( J- ) = ( - 1 ) (x wyil + XllVl0 ) + . 2 ( _^ iq ) (2.7) 

The solvability condition requires that the e lw * part of the right hand side be 
orthogonal to the left eigenvector of L. This leads to 

- 2 = (2-8) 
The perturbative result for uj up to 0(A 2 ), after setting A = 1 is 

- = 1-^ (2-9) 

where A is the amplitude of the limit cycle. If xq and yo be the initial values of 
x and y, then we can write 

w = 1 _( X0 -l? + (V0-lf 

12 v ; 

We have checked this result numerically. The result are shown in Fig pp. The 
good agreement between the computed frequencies and the obtained from Ea H2.10l) 
is apparent. 



3 



0.982 



0.985 



0.988 



0.991 



0.997 



0.994 s 



1.00 




1.003 



o:b Eras Era um 



TTOS T7T 



T35 T.2 



Figure 1: Initial number of population density modifies its periodicity. Initial 
potulation density (x=y) is plotted with corresponding frequency (lo) of its 
oscillation. 

3 The modulated system at low and moderate 
frequencies 

With the help of the amplitude dependent frequency of the previous section, 
we explored Lotka-Volterra population dynamics under external drive. These 
systems have been extensively studied under effect of random perturbation 
whereas under periodic forcing they were studied much less ^7] . We explored it 
in presence of periodic perturbation. Periodic perturbation leads to diurnal or 
annual influences on the predator-prey systems and leads to periodic intrinsic 
growth rate in the prey population. 
Under periodic forcing, the system [Ea i|2.1[l ] becomes 



We note that if e = then the solutions can be written as periodic trajectories 
around (1,1) with a frequency that is dependent on the initial conditions. For 
e = we have equivalent oscillator 



x — x — xy + ex cos Qt 

y = -y + xy 



(3.1) 



Shifting to the (x\,y\) variables, we have 



x\ = —yi — xiyi + e(l + X\) cos fit 
y\ =x x + x x yi 



(3.2) 



.Ei = -uyi 



(3.3) 



2/i = wxi 
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with u> given by Ea H2.1U|) . In this technique of equivalent linearisation Ea l|3.3[) 
can be rewritten as; 



x\ = —ujyi + e(l + xi) cos fit 
y'l = LUXl 



(3.4) 



Eliminating y\ we have 

x\ — ecos(ilt)x\ + (uj 2 + eft sin Qt)xi — — eilsinSlt (3-5) 
Defining Xl = z - then, 

, + ^-e-(,cosOi-^-^ y ) = (3.6) 

When the driving amplitude is small, (e 2 <C e) 

I - ecosOii + (w 2 + e£lsinf2<)z = (3.7) 

Using the transformation z = eafr slnSlt p(i), we arrive at 

p + (lu 2 + en sin Qt)p = (3.8) 

correct to 0(e). We recognize Ea (|3.8Jl as a Mathieu equation and note that 
there will be a periodic response at a frequency of § provided we fulfill the 
condition 

w = ~(l±e) or 

(xo-i) 2 + ( yo -i) 2 n (3 ' 9) 

12 2 1 J 

This is in effect a construct on the initial conditions. The result found in Ea l|3.9JI 
implies that for every modulating frequency f2, there will be some initial condi- 
tions for which a periodic response will be possible. For initial conditions in the 
range §(l-e) < 1- (. x o-i) 2 Myo-i) 2 < 0(i + e ), the response will be unbounded 
while for initial conditions outside this range the response will be bounded and 
in general quasi-periodic - the periodic motion resulting when Ea (|3.9l) is satis- 
fied. 

In Fig|3] , the unboundcdncss of the solution is shown for fl = 2. The width of 
the unbounded region is exactly 2e. The original set of equation have nonlin- 
ear terms that prevents the dynamics becoming unbounded. Nevertheless we 
get the signature of unboundedness through the steep increase of the width (77) 
FigjBl of the phase space trajectory in the same frequency range. r\ is poten- 
tially very different from the width £ we discussed later. 77 is a effect of the 
dynamics of two competing frequencies having resonating effect at comparable 
values, whereas £ shows the presence of two different order of timescale. 
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Figure 2: Time series (a) of prey population (xi) shows the signature of un- 
boundedness of the linearized solution. Fig© Change of width (77) (a,b) reflects 
the resonance in linearized rezime. As we move from linearized unbounded so- 
lutions 77 decreases (c,d). solid lines (a,c) represent the unperturbed trajctories. 
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Figure 3: a) Bounded and unbounded solution from Mathieu equation [Eq(??)] 
campares well with the signature of resonance in original dynamics (b). In (b) 
77 is plotted with initial conditions of population density (xq = yo). 



4 High Frequency Limit 

In this section, we consider the variation in the parameter of the Lotka Volterra 
model to be very rapid, i.e. the period Q of the forcing in Ea l|4.1|l . is much 
greater than the frequency of the unforced system. The system fluctuates 
rapidly from the unperturbed trajectory and the phase space trajectory is broad- 
ened Fig[2]. This broadening (£) depends mainly on the frequency of forcing 
and feebly on forcing amplitude. As O increases this broadening decreases very 
rapidly. The analysis of the system follows a procedure explained by Landau 
and lifshitz. We split the variables X\ and y\ into two parts, 

*i = Xi + 6 
yi = Yi + & 

where £1, £2 carry the rapid variations (scale of O) and the averaged quanti- 
ties X\ and Y\ carry the slow variation (order unity) of the original model. 
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Figure 4: (a,c)Time series shows how the amplitude (£1) of high frequency 
perturbation changes as we vary frequency of perturbation (f2) from 10 (a) to 
30 (b). Phase space plots (b,d) also show these through its change of width, 
the solid lines represent the unperturbed trajectory. 



Substituting in Ea (|4.1jl : 

Xt+fr = -Y 1 -t2-(X 1 +Z 1 )(Y 1 +Z 2 ) + e(l + X 1 +Z 1 )caant (4.2) 

= + + ( 4 - 3 ) 

In this section, we do not restrict e to be small. Instead, we note that since 
£ ~ cos fit, typically £ ^> £. We choose the dynamics of £1,2 to be 

= eX 1 cosfit- X^ 2 (4.4) 
6 = Yifr (4.5) 

In the above X x and Yj. may be treated as constants on the scale of variation 
over a period Q,^ 1 . In that situation, Eas (|4.4l) and l|4.5l) reduce to 

6 + YxXx^x = -eX x sin fit (4.6) 

Since f2 S> \fX\Yx, we write the approximate solution of Ea l|4.6|l as 

ei = ^sinm (4.7) 

We now average Eas <|4.2|) and i|4.3[l over the fast variations, drop terms that are 
0(fl^ 2 ) and obtain 

*\—*-™ ,4.8, 

r, = x 1 + x i k, 

which is the unforced system in the coarse-grained variables. The picture, which 
emerges is as follows: under a high frequency modulation, the motion can be 
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explained as a rapidly oscillatory motion which is virtually identical to the 
trajectory of the unforced system. The amplitude of the fast variations decays 
according to ^ as shown in Ea (|4.7|) . Numerical simulation of the system bears 
out these expectations. Fig[S] 




Figure 5: width £1 is plotted against the frequency of modulation (SI). The 
fitting function (F(x)) supports the form of analytical result. 



5 conclusion 

In this paper we have shown that under lotka-volterra dynamics populations of 
predator and prey depend on their initial population densities. This dependence 
of initial population brings forth changes in dynamics but helps us formulate 
an effective Mathew-Hill type equation for this 2-species system under periodic 
seasonal changes. This also indicates slow seasonal changes have more effects in 
the variation of population than its faster counterpart. 

References 

[1] A.J Lotka, J. Am. Chem Soc. 42, 1595 (1920). 

[2] V. Voltcrra, Atti R. Accad, Naz. Lincei, Mem. CI. Sci. Fis,. Mat. Nat. 2, 
31 (1926). 

[3] H. Haken, Synergetics, 3rd ed. (Springer- Verlag, New York, 1983). 

[4] D. Neal, Introduction to Population Biology (Cambridge University 
Press,Cambridge, U.K. 2004). 

[5] R. M. May, Stability and Complexity in Model Ecosystems (Princeton Uni- 
versity Press, Princeton, 1973). 



8 



[6] J. Maynard Smith, Models in Ecology (Cambridge University 
Press, Cambridge, U.K. 1974). 

[7] J. D. Murray, mathematical Biology, 3rd ed. (Springer- Verlag, New York, 
2002), Vols. I and II. 

[8] T. Fukai and S. Tanaka, Neural Comput. 9, 77 (1997). 

[9] Zhang Yi,K. K. Tan, Phys. Rev. E, 66, 011910 (2002). 

[10] Luca Bianco, Federico Fontana, Vincenzo Manca 

[11] J.J. Hopeld, Proc. Natl. Acad. Sci. USA 81, 3088 (1984). 

[12] K. Gopalsamy and X.Z. He, Physica D 76, 344 (1994). 

[13] M. Mobilia, T. Goeorgiev, U. C. Tauber, Phys. Rev. E, 73, 040903(R) 
(2006) 

[14] A. R. Balsara, K. Lindcnberg and K. E. Shulcr, J. Stat. Phys,27 787 (1982) 
[15] A. B. Budgor, B. J. West Phys. Rev. A, 17 370 (1978) 
[16] Mikhail F. Dimcntberg, phys. Rev. E, 65 036204 (2001) 
[17] Sabin G.C.W, Summer D, Math Biosci. 113 91 (1993) 



9 



